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We recently proposed a new structural model for the Si(331)-(12x 1) surface reconstruction con- 
taining silicon pentamers as elementary structural building blocks. Using first-principles density- 
functional theory we here investigate the stability of a variety of adatom configurations and deter- 
mine the lowest energy configuration. We also present a detailed comparison of the energetics be- 
tween our model for Si(331)-(12x 1) and the adatom-tetramer-interstitial model for Si(110)-(16x2), 
which shares the same structural building blocks. 
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INTRODUCTION 



High-index silicon surfaces are of interest not only 
from a fundamental point of view, but also as poten- 
tial substrates for electronic device fabrication. The 
reconstructed Si(331)-(12xl) surface is of particular 
importance since it is the only confirmed planar silicon 
surface with a stable reconstruction located between the 
(111) and (110) directions. Due to its large unit cell 
and its pronounced structural anisotropy it represents 
a promising template for the growth of low-dimensional 
heteroepitaxial nanostructures such as quantum dots 
and nanowiresii 

Since its discovery more than 17 years ago^ several 
structural models containing dimers and adatoms 
as elementary structural building blocks have been 
proposed. '^''^ We recently revealed the presence of an 
additional building block on the Si(331)-(12xl) surface?^ 
Using scanning tunneling microscopy (STM) we were 
able to resolve for the first time rows of pentagon 
pairs running across the Si(331) surface. Very similar 
pentagons were observed earlier on the reconstructed 
Si(llO) and Ge(llO) surfacesi^i^ Inspired by struc- 
tural elements encountered on Si(113) and Ge(113) 
surfaces r''^'^ An et al^ have proposed an adatom- 
tetramer-interstitial (ATI) model for the Si(110)-(16x2) 
reconstruction containing silicon pentamers as building 
blocks which explain the pentagons observed in STM 
images. The stability of the ATI model has subsequently 
been tested theoretically by means of first-principles 
total energy calculations .-^Ji Based on a detailed 
analysis of our experimental results and a comparison 
between the Si(llO) and Si(331) surface, we proposed a 
new structural model for the Si(331)-(12x 1) reconstruc- 
tion containing silicon pentamers as essential building 
blocks.^ 

In order to account for the pentagons observed in 
our STM images, the model contains two pentamers 
per (12x1) unit cell, which saturate a certain number 
of energetically unfavorable dangling bonds of the 
bulk-truncated Si(331) surface. Some of the remaining 
dangling bonds are saturated by simple adatoms. The 



pentamer positions and bonding configurations were 
determined exclusively using experimental information 
and by comparing Si(331) to Si(llO) as discussed 
in detail in Ref. [g. Although the available body of 
experimental results allows to narrow down the location 
of the adatoms, currently we do not have sufficient 
evidence to precisely determine the adatom positions 
on the surface solely based on information derived from 
experiment. In this paper we use first-principles total 
energy calculations to determine the precise adatom 
positions by finding the lowest energy configuration 
among various alternative adatom configurations. 
The paper is organized as follows: After briefly review- 
ing the methods used for the calculations, we focus 
on the bulk-truncated Si (331) surface and compare its 
energetics to other silicon surfaces, which are known to 
be stabilized by surface reconstructions. In a second step 
pentamers are introduced and the consequences on the 
total energy are discussed. Then we consider a variety of 
candidate adatom configurations and determine the one 
with the lowest energy. At the end we discuss the im- 
pact of the remaining dangling bonds on the total energy. 



II. THEORY 
A. Computational details 

Our results are based on first-principles calculations 
within density- functional theory (DFT) )^^i^^ using the 
local-density approximation (LDA) for the exchange- 
correlation functionalJ^ii^ We adopt a plane-wave ba- 
sis set and norm-conserving pseudopotentials, as im- 
plemented in the PWSCF code of the QUANTUM- 
ESPRESSO suite. ^° Pseudopotentials for sihcon and hy- 
drogen are choosen in the von Barth style. ^^'^^ 
The Si (331) surface was simulated by means of an atomic 
slab and a repeated supercell, in order to recover a three- 
dimensional periodicity. The slab thickness corresponds 
to 10 sihcon double layers (DLs), and the periodic images 
are separated by 9.4 A of vacuum. Since each bulk DL 
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TABLE I: Surface energies for unrelaxed bulk-truncated surfaces of various orientations. The numbers in parenthesis give the 
number of surface atoms and dangling bonds per (1x1) unit cell respectively. 



Surface 
orientation 


Surface atom 
density (A"^) 


Dangling bond 
density (A"^) 


Estimated surface 
energy (meV/A^) 


Calculated surface 
energy (meV/A^) 


(111) 


0.078 (1) 


0.078 (1) 


90.1 


113.6- 


(100) 


0.068 (1) 


0.136 (2) 


155.9 


149.2*^ 


(113) 


0.082 (2) 


0.123 (3) 


141.0 


137.9i^ 


(110) 


0.096 (2) 


0.096 (2) 


110.3 


127.3^ 


(331) 


0.093 (3) 


0.093 (3) 


107.3 


127.2 " 
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contains 24 silicon atoms, our supcrccU counts 240 silicon 
atoms plus a number of silicon atoms at the top surface 
(varying from to 14 depending on the specific recon- 
struction) . The bottom surface of the slab is hydrogen- 
passivated (with 36 hydrogen atoms). 
Since we are interested in small energy differences, we 
adopt strict convergence criteria with respect to both the 
kinetic-energy cutoff and the sampling of the surface Bril- 
louin zone. In particular, the inclusion of plane waves up 
to 544 eV (40 Ry) and the use of a xy shifted 3x4x1 
Monkhorst-Pack grid^^ for the k point sampling allowed 
us to ideally compute total energy differences with an ac- 
curacy of the order of 15 meV/supercell. By taking into 
account the fact that our relaxed structures still show 
small residual forces (see below), we estimate our numer- 
ical accuracy on calculated surface energies to be better 
than 0.1 meV/A^. All structures were optimized till the 
largest residual force on the mobile atoms was less than 
0.026 eV/A (0.001 Ry/Bohr). 

The hydrogen atom positions were determined in a first 
run by optimizing them on a bulk-terminated slab with 
fixed bulk-like positions for silicon atoms, with top and 
bottom surfaces saturated by hydrogen. In all the follow- 
ing calculations we have taken hydrogen atomic positions 
at the bottom surface as fixed. In a second step, we in- 
spected the decay of forces in the direction perpendicular 
to the surface, for a slab containing a non-relaxed surface 
reconstruction. A nearly force-free bulk region is reached 
after the five topmost DLs. Thus, only the topmost five 
bulk silicon DL were allowed to relax together with the 
surface (reconstruction) atoms, where the inner DLs have 
been kept fixed. All results presented below are based on 
a slab relaxed in this way. 



B. Determination of the surface energy 

The surface energy i?surf of a double-sided symmetric 
slab containing iVsi silicon atoms can be written asi^ 



surf = -(£^tot(^Si) -MSiA^Si) 



(1) 



value can be obtained from a simple bulk silicon calcu- 
lation. -Etot(A^Si) is the total energy of the slab. The 
factor of 2 comes from the fact that there are two iden- 
tical surfaces. For a non-symmetric slab for which only 
the bottom surface is saturated by hydrogen atoms, we 
have 

£^surf + £^s^.rf - ^^tot [Nsi, Nu) ~ MSi^^Si ^ /^hA^H (2) 

Here, E^^.f is the surface energy of the hydrogen satu- 
rated surface, /zh is the chemical potential for the hydro- 
gen atoms, and i?tot(A'^si, A'^h) is the total energy of the 
slab whose bottom surface is saturated by A^h hydrogen 
atoms. Therefore, we have 



£^surf = £^tot(A^si, A^h) - MsiA^si - mhA^h - El 



(3) 



In order to avoid the determination of /in, we consider 
the following: The total energy of the slab can be written 
as 



EtotiNsuNn) = ^Efot + ^Et, 



(4) 



where /isi is the chemical potential of silicon, being sim- 
ply the bulk energy per atom at zero temperature. This 



where E^^. is the total energy of a symmetric hydrogen 
saturated slab and Etot is the total energy of a symmetric 
slab containing the surface reconstruction. The value 
E^^^ has been calculated in the preliminary run used for 
the determination of the hydrogen positions. The surface 
energy (per unit area) for the top surface carrying the 
surface reconstruction can then be determined using 

7 = Esuri/A = (£;tot(A^Si, A^h) - ^E^ot - Nsi^isi)/A (5) 

where N^i is the number of silicon atoms in the upper 
half slab and A the surface area of the slab. 



III. RESULTS AND DISCUSSION 

A. Bulk-truncated surface 

In order to define an energy reference for the follow- 
ing discussion we have determined the surface energy for 
the unrelaxed bulk-truncated Si (331) surface shown in 
Fig. [B^a) and (b). The value of 127.2 meV/A^ may be 
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compared in Tab. U with the values for other surface ori- 
entations. 

The surface energy of a bulk-truncated surface can be 
estimated from the cohesive energy of silicon, which is 
Ec — 4.6 eV per atom.^'* The cohesive energy of a solid 
is the energy per atom required to break the atoms of 
the solid into isolated atomic species. Since a bond is 
formed by two atoms and each silicon atom makes four 
bonds, the energy required to break a bond is half the 
cohesive energy, i.e. Ef, = 2.3 eV per bond.^^ The surface 
energy 7 can now be estimated by multiplying the bond 
energy Ei, by the density of broken bonds n.^^- The latter 
is the number of dangling bonds N per surface area A, 
n = N/A. Thus 7 = NEt/2A = nEb/2. The factor of 2 
takes into account that two surfaces are created. 
Note that the dangling bond density n is not necessarily 
equal to the surface atom density (see values in paren- 
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thesis in Tab. This is due to the fact that a surface 
atom may either carry one or two dangling bonds. (111)- 
like surface atoms carry a single dangling bond, whereas 
(lOO)-like surface atoms carry two dangling bonds (see 
Fig. 3(a) and 2(a) respectively in Ref. 27). On the (113) 
surface, both types of surface atoms coexists, whereas on 
the (110) and the (331) surface only (lll)-type surface 
atoms occur (see Fig. 5(a), 6(a) and 7(a) respectively in 
Ref. [13). 

In Tab. U we summarize estimated and calculated surface 
energies of unrelaxed bulk-truncated surfaces for a series 
of surface orientations. Although the simple estimate 
based on the cohesive energy fails to precisely reproduce 
the values determined with the help of first-principles 
methods, it correctly reproduces the energy ordering of 
the various surfaces. The estimated values tend to under- 
estimate the surface energy for surfaces containing only 
(lll)-type surface atoms, whereas they overestimate the 
surface energy for surfaces with (lOO)-type surface atoms. 
For Si(113), which contains both types of surface atoms, 
these two tendencies compensate resulting in the best 
agreement between estimated and calculated value. 
From the theoretical point of view, the most stable sur- 
face is the one with the lowest surface energy. Conse- 
quently only structural models for the reconstructed sur- 
face with a surface energy lower than the bulk-truncated 
surface must be considered. Allowing the bulk-truncated 
surface to relax and determining the corresponding sur- 
face energy, we obtain an even more stringent condition 
on the surface energy for any valid candidate model. The 
relaxed structure of the bulk-truncated Si(331) surface is 
shown in Fig.[2](a). Note that although the relaxation was 
performed within the (12x1) supercell, it did not break 
the (1x1) periodicity within our numerical accuracy. The 
surface energy for the relaxed bulk-truncated surface is 
given in Tab. |TT] along with the surface energies for other 
structures for the Si(331) and Si(llO) surface discussed 
hereafter. Relaxation of the bulk-truncated Si(331) sur- 
face within the (12x1) unit cell results in a reduction of 
the surface energy of 16.3 meV/A^, a value to be com- 
pared with the corresponding value of 21.2 meV/A^ for 
Si(llO) obtained by Stekolnikov et alt^ via relaxation 
of the bulk-truncated Si(llO) surface within the (16x2) 
cell. 



FIG. 1: (Color online) Side (a)/(c) and top view (b)/(d) 
of the unrelaxed bulk-truncated Si(331)-(12x 1) surface unit 
cell/unrelaxed 'pentamers only' model. The bulk directions 
are given. In (b) / (d) only the topmost atoms contained within 
the rectangle in (a) /(c) are shown. Subsurface atoms are 
shown in grey, undercoordinated surface atoms in black, pen- 
tamer atoms in red. In (a) the Si(lll) and Si(llO) terraces 
composing the Si(331) surface are indicated. In (b) the (1x1) 
and (12x1) surface unit cell vectors ai, and Ai, A2 re- 
spectively are also drawn. The orientation of the glide plane 
is indicated by a dashed line. In (b) the four different types 
of adatom adsorption sites T4,, H3, A3 and A'-^ are represented 
by markers. In (d) the 2x12 possible adatom positions after 
introduction of the two pentamers are shown. 



B. Pentamers 

The bulk-truncated Si (331) surface consists of alter- 
nating single and double rows of undercoordinated sur- 
face atoms running along the [110] direction (black atoms 
in Fig. [Ija) and (b)). Each of the 36 (lll)-type surface 
atoms per (12x1) unit cell carries one dangling bond, 
which renders the unreconstructed structure energeti- 
cally unfavorable. 

In our model, 10 out of the 36 dangling bonds are satu- 
rated by two silicon pentamers observed as pentagons in 
STM images. Their integration onto the bulk-truncated 



4 



a) bulk-truncated b) pentamers only 
relaxed Y=11 0.9 y=114.9 




c)T4 7=111.8 d)H3Y=112.8 




e)A3Y=112.7 f)T4ay=111.2 




FIG. 2: (Color online) (a) Top view of the relaxed bulk- 
truncated Si(331)-(12x 1) surface, (b) Top view of the relaxed 
'pentamers only' model, (c)-(j) Various candidate adatom 
configurations of the adatom-pentamer model. The surface 
energy 7 for each model is given in units of me 



Si (331) surface was discussed in detail in Ref. It is 
important to note that the bulk-truncated Si(331) sur- 
face may be viewed as a highly stepped surface consist- 
ing of small alternating Si(lll) and Si(llO) terraces (see 



TABLE IL Comparison between the surface energies of vari- 
ous surface structures of Si(331) and Si(llO) 





Surface 


energy (moV/A^) 


Structure 


Si(331) 


Si(llO) 


bulk-truncated 


127.2" 


127.3i^ 


bulk-truncated relaxed 


110.9" 


lOe.li^iii 


pentamers only 


114.9" 


los.eHiii 


adatoms and pentamers 


107.9" 


103.SHiii 
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Fig. IHa)). This allows to anchor the pentamers on the 
Si(331) surface in exactly the same local binding config- 
uration as on the Si(llO) surface, where pentamers are 
also observed. Introducing two pentamers per unit cell 
as shown in Fig. [ijc) and (d) reduces the number of dan- 
gling bonds from 36 to 26 and lowers the surface energy 
(for the corresponding relaxed model shown in Fig.j^Jb)) 
by 12.3 meV/A^ with respect to the unrelaxed bulk- 
truncated surface (see Tab.ln|. This value is comparable 
with the reduction of 18.7 meV/A^ obtained when intro- 
ducing four pentamers onto the Si(110)-(16x2) surface 
(corresponding to the tetramer-interstitial (TI) model in 
Ref. JJ,, note that in the TI model the step is already 
included into the unit cell accounting for 1.0 me en- 
ergy reduction). However, the 'pentamers only' models 
for both, Si(331) and Si(llO), are not stable structures, 
since the relaxed bulk-truncated surface for both surfaces 
possesses a lower surface energy. This is not surprising, 
since the introduction of pentamers suggested by STM 
images, leaves a lot of dangling bonds unsaturated. So 
an additional building block is required to stabilize the 
structures containing pentamers. 

C. Adatoms 

In the ATI model for Si(110)-(16x2), the pentamer 
structure is stabilized by adatoms. Among the various 
elementary structural building blocks observed in sili- 
con surface reconstructions (see Ref. i27, for a review), 
adatoms impose themselves also as the best choice for 
Si (331) to saturate some of the remaining 26 dangling 
bonds for the following two reasons: Firstly, as men- 
tioned before, the bulk-truncated Si(331) surface may be 
viewed as being composed of small alternating (110) and 
(111) terraces. Since adatoms are observed on both the 
Si(llO) and the Si(lll) surface, the Si(331) consequently 
provides exactly the same local environment for adatoms 
as on Si(llO) and Si(lll). Secondly adatoms allow to 
explain the additional protrusions seen in STM images, 
which can not be attributed to the pentamer atoms?^ 
An adatom requires three unsaturated surface atoms to 
be attached to. On the (111) terraces of the Si(331) sur- 
face, adatoms may occupy two possible sites indicated in 
Fig. mb) . These geometries are distinguished as hollow 
(i/a) and atop {T4) sites depending on whether the sub- 
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strate atoms below the adatom is found in the fourth or 
second layer. In sites the adatom is three-fold coordi- 
nated, in T4 sites the adatom is approximately four-fold 
coordinated due to the substrate atom directly below in 
the second layer. 

On the (110) terraces of the Si(331) surface, we can also 
distinguish between two different three-fold coordinated 
adsorption sites depending on their orientation labeled 
A3 and ^3 in Fig.l^b). Note that on the bulk-truncated 
(110) surface these two sites become equivalent due to 
the mirror symmetry of the (110) surface along the [110] 
direction (not shown). 

Careful inspection of Fig.fTfd) shows that after introduc- 
ing two pentamers per unit cell there are exactly 2x12 
possible positions for placing adatoms, the factor of 2 is 
due to the glide plane symmetry along the [116] direction 
(dashed line in Fig.[ljd)). Out of the 2x12 candidate ad- 
sorption sites, we identify 2x3 T4 sites, 2x4 H3 sites, 2x3 
A3 sites and 2x2 A3 sites. This however does not imply 
that a maximum of 24 adatoms can be placed on the sur- 
face. The reason for this is as follows. After placing a 
first adatom on the surface, some surrounding candidate 
sites can not be occupied by an adatom anymore, since 
they shared the same dangling bond with the now occu- 
pied site, whose adatom saturates this dangling bond. 
A first pair of adatoms must be placed between two suc- 
cessive pentamers, in order to account for the clear pro- 
trusions seen in STM images.^ Due to the glide plane 
symmetry adatoms must always occur in pairs. There 
are only three candidate sites for this first adatom: a 
T4 site (Fig. M.C)), a H3 site (Fig. ^d)) and a A3 site 
(Fig. mje)). Using first-principles calculations we have 
tested these three possible adatom configurations and 
found that the T4 site is the most stable of the three 
(surface energies for each model are also given in Fig. [2]). 
This is not surprising, since T4 sites are commonly fa- 
vored over the H3 sites on (111) surfaces such as in the 
Si(lll)-(7x7) reconstruction and the Gc(lll)-c(2x8) re- 
construction (see Ref. [13 for a review of their structural 
models). In addition the T4 adatom position reproduces 
best the protrusion observed by STM between two suc- 
cessive pentagonsi^ 

Consequently, in order to reduce the large number of 
potential candidate adatom configurations to a compu- 
tationally manageable number of models, we have re- 
stricted our further investigations to the subset of models 
containing the T4 adatom pair. In a second step we have 
investigated all possible remaining positions for a second 
pair of adatoms. There are exactly five candidate con- 
figurations to be considered (shown in Fig. [21^f)-(j)). Us- 
ing first-principles calculations we find that a A3 adatom 
position (configuration T4e) gives the most stable con- 
figuration. The position of the A3 adatom pair is also in 
agreement with light protrusions seen in STM images^. 
After placing this A3 adatom pair there is no more can- 
didate site left on which one could place an additional 
adatom. 

As discussed before, the introduction of two pentamers 



TABLE III: Surface energies for silicon surface reconstruc- 
tions 



Reconstruction 


Surface energy 
(meV/A^) 


Si(lll)-(7x7) 


84.9''* 


Si(100)-c(4x2) 


SS.O'"* 


Si(113)-(3x2) 


87.92V87.4ii 


Si(110)-(16x2) 


lOS.SHiii 


Si(331)-(12xl) 


106.4" 
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is not sufficient to obtain a lower surface energy than for 
the relaxed bulk-truncated surface. Inspection of Tab. HIl 
shows that adding adatoms to the pentamers allows to 
sufficiently lower the surface energy to render the re- 
constructed surface more stable than the relaxed bulk- 
truncated surface, for both the Si(331) and the Si(llG) 
surfaces (reduction of surface energy by 7.0 meV/A^ and 
4.8 meV/A^ respectively with respect to the 'pentamers 
only' models). 



D. Rest atoms 

Introducing two T4 and two A3 adatoms per unit 
cell into our structural model for the Si(331)-(12x 1) 
reconstruction further reduces the number of dangling 
bonds from 26 to 14 (10 of the 36 dangling bonds are 
saturated by the pentamers, introducing 4 adatoms 
further eliminates 4x3 dangling bonds (here we do not 
count the dangling bond carried by each adatom)). The 
remaining 14 surface atoms, which still carry a dangling 
bond, are so called rest atoms (black atoms in Fig. [Hi)). 
Rest atoms are common structural building blocks 
accompanying adatoms. The adatom model for the 
Ge(lll)-c(2x8) reconstruction is built from four 
adatoms and four rest atoms. For this case, we have an 
ideal one-to-one correspondence between adatoms and 
rest atoms. Each adatom delivers one electron from its 
dangling bond to fill the dangling bond of the rest atom 
resulting in a further decrease of the surface energy. 
In contrast, in the famous dimer-adatom-stacking fault 
(DAS) model of the Si(lll)-(7x7) reconstruction, there 
are 12 adatoms, but only 6 rest atoms plus the corner 
hole atom, which behaves similarly to a rest atom. 
In this case, using a simple ionic picture, we have 7 
electrons which are transferred from the adatoms into 
the rest atom plus corner hole atom states, but five 
electrons remain in the adatoms bands. For the ATI 
model of the Si(110)-(16x2) reconstruction, there are 8 
adatoms and 12 rest atoms. So even when transferring 
the 8 electrons from the adatoms to the rest atoms, only 
a fraction of the rest atoms may saturate their dangling 
bonds. Similarly in our model, with four adatoms and 
14 rest atoms, we expect some remaining unsaturated 
rest atoms. 
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The unequal adatomirest atom ratio between the ATI 
model (8:12=0.67) of Si(llO) and our model (4:14=0.29) 
of Si(331) might on the first view explain why the ATI 
model exhibits a slightly lower surface energy than our 
model (energy difference 4.1 meV/A^, see Tab. 
Reinspection of Tab. HIl however reveals that this energy 
difference between Si(331) and Si(llO), whose corre- 
sponding unrelaxed bulk-truncated surfaces have almost 
the same energy (energy difference only 0.1 meV/A^), 
is already present for the respective 'pentamers only' 
models (energy difference between 'pentamers only' 
models 6.3 meV/A^). It is also interesting to note 
that the introduction of the adatoms leads to a larger 
reduction of the surface energy for Si(331) (7.0 meV/A^) 
than for Si(llO) (4.8 meV/A^) with respect to the 
'pentamers only' models. 

Introducing hypothetically an additional adatom pair 
into our model for Si(331) would improve its adatom:rest 
atom ratio from 4:14=0.29 to 6:8=0.75. The introduc- 
tion of the first pair of adatoms (T4) accounts for 
a reduction of 3.1 meV/A^ of the surface energy. 
Adding a second pair of adatoms {A3) to the first pair 
allowed to reduce the surface energy by 3.9 meV/A^. 
Following this tendency, by assuming an energy gain 
of 3-4 meV/A^ for the introduction of a hypothetical 
third pair of adatoms into our model, suggests that the 
surface energy of Si(331) comes to lie close to the value 
determined for Si(llO) by Stekolnikov et al.^^. However, 
note that the rest atom geometry of T4e does not allow 
the introduction of an additional adatom pair. When 
considering the subset of T4x models (with x=a-e), 
only configurations T4a and T4b admit a pair of free 
adatom sites, resulting in configuration T4a-|-b. But this 
configuration was discarded because it is not compatible 
with the experimental STM images.- 
Several more three adatom pair configurations are 
possible when studying II3x (8/7 possible two/three 
adatom pair configurations) or A3x configurations (5/1 
possible two/three adatom pair configurations). Even 
more possibilities may exist when allowing the placement 
of the first adatom pair on a different site. Here we 
focused on the family of T4x models, whose selection was 
motivated and justified in section C. Before testing other 
potential configurations, further experimental hints 
on the position of the adatoms from complementary 
experimental techniques are necessary to assist in the 
search for the lowest energy configuration. 
When compared to the surface energies of reconstruc- 
tions on surfaces with other orientations listed in 
Tab. nH the models for the Si(110)-(16x2) and the 
Si(331)-(12xl) reconstruction have surface energies 
which are of the order of 20 meV/A^ higher than 
the surface energies for Si(lll)-(7x7), Si(100)-c(4x2) 
and Si(113)-(3x2). As it was already mentioned by 
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Stekolnikov et al}^ in the context of the ATI model 
for Si(110)-(16x2), a completely new, more complex 
reconstruction model cannot be excluded. Nevertheless, 
the advantage of our method for developing a structural 
model for Si(331)-(12x 1) consists here in integrating 
elementary structural building blocks such as pentamers 
and adatoms identified on the (110) and (111) silicon 
surfaces^^ onto the bulk-truncated Si(331) surface. The 
consistency of the pentamer building block with all 
currently available experimental data on Si(331)-(12x 1), 
but also on Si(110)-(16x2), represents further indirect 
evidence in favor of the two structural models. It should 
be noted, that our model for Si(331) is far from a simple 
extension of the Si(llO) surface to its vicinals. The 
pentagons on Si(llO) (see Fig. 6 in Ref. '27') are arranged 
in pairs pointing in opposite directions, whereas on 
Si (331) pairs point in the same direction. Also the 
direction of the double pentagon rows appears to be 
completely unrelated between the two surfaces. Using 
pentamers and adatoms we succeeded in combining 
these seemingly confiicting experimental observations in 
a coherent picture. 



IV. SUMMARY AND CONCLUSIONS 

Using first-principles density-functional theory we 
studied the stability of the structural model for the 
Si(331)-(12x 1) surface reconstruction, which we recently 
proposed. In particular we determined the energetically 
most favorable adatom configuration among a variety of 
candidate structures. We find that a configuration with 
a pair of T4 adatoms and a second pair of A3 adatoms 
results in the lowest energy. We also discussed the im- 
plications of the rest atoms on the energy and compared 
our model to other silicon surface reconstructions. 
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